Magnetic excitations in the low-temperature ferroelectric phase of multiferroic 
YMn 2 5 using inelastic neutron scattering 
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We studied magnetic excitations in a low-temperature ferroelectric phase of the multiferroic 
YMn205 using inelastic neutron scattering (INS). We identify low-energy magnon modes and es- 
tablish a correspondence between the magnon peaks observed by INS and electromagnon peaks 
observed in optical absorption [l[. Furthermore, we explain the microscopic mechanism, which re- 
sults in the lowest-energy electromagnon peak, by comparing the inelastic neutron spectral weight 
with the polarization in the commensurate ferroelectric phase. 

PACS numbers: 78.70.Nx, 75.30.Ds, 75.85.+t, 78.20.-e 



The coupling of spins with electric polarization in mul- 
tiferroics gives rise to a wide variety of magnetoelectric 
effects actively studied over the past ten years In 
many multiferroics, electric polarization is induced by 
magnetic structures formed when conventional magnetic 
orders are frustrated by lattice geometry and competi- 
tion between spin interactions. Magnetic frustration can 
yield unusual excitations, such as spinons, monopoles and 
localized non-dispersive modes [7HlQj. These magnetic 
excitations can couple to an electric field, which results 
in the electromagnon peaks in optical absorption [l], [TlJ • 

Many spectacular magnetoelectric phenomena have 
been found in manganese oxides, such as mixing of 
magnons with acoustic phonons, clamping between fer- 
roelectric and magnetic domain walls, and unusual 
magnetoelectric vortices in hexagonal YMn03 [lil-lli). 
In orthorhombic 74Mn03 compounds with ^4=Tb, Dy, 
(Eu, Y), an incommensurate spiral spin order induces 
electric polarization through the relativistic 'inverse 
Dzyaloshinskii-Moriya mechanism' 15- These com- 



pounds show magnetically-induced polarization flops, gi- 
ant magnetocapacitance and electromagnon peaks in op- 
tical absorption [H, Qj], Qjl, [lo|. Smaller cations on A- 
sites, e.g., Ho, favor a commensurate collinear magnetic 
order whose polarization, induced by the non-relativistic 
Heisenberg magnetoelectric coupling, is an order of mag- 
nitude larger than that of the spiral state [21] . 

A similar competition between collinear and non- 
collinear ferroelectric states is found in orthorhombic 
^M^Os compounds. In this family, electric polarization 
is also larger in collinear states, while electromagnons and 
giant magnetocapacitance are observed in non-collinear 
ones [3, 13[ UK • In YMn 2 5 , Mn spins order at 44 K in 
a periodically modulated structure with an incommen- 



surate wave vector (q x ,0,Qz)- This state is paraelectric. 
Upon further cooling by a few degrees, first q z and then 
q x locks to a commensurate value and this commensu- 
rate (CM) phase, having the wave vector qc = (f ,0, |), 



is ferroelectric [23|, [24[. At 19 K the wave vector of 



the magnetic structure becomes again incommensurate: 
qic = (0.48,0,0.288). This low-temperature incommen- 
surate (LTI) state is also ferroelectric but has a much 
smaller polarization 23], [25| . This unusual sequence of 
transitions results from geometrical frustration in the Mn 
spin lattice [HUHI. 



Magnetic orders in the CM and LTI phases were for 
some time a matter of controversy. However, several 
groups recently reported similar spin structures [H, [30l — 
132] . In the CM phase spins are approximately collinear 
and lie nearly antiparallel along the a axis, due to the rel- 
atively strong antiferromagnetic (AFM) exchange along 
this direction. The AFM chains are also discernable in 
the LTI state, which, however, is strongly non-collinear 
as the angles between spins in neighboring AFM chains 
are close to ±90°. The LTI spin structure, which can 
roughly be characterized as a spiral where spins rotate in 
the ab plane, is actually rather complex as it also contains 
considerable be and ac spiral components 23|, [3l|, l32| . 



A phenomenological analysis of the observed spin orders 
was used to show that the spontaneous electric polariza- 
tion in the CM phase largely results from non-relativistic 
exchange interactions, while the much smaller polariza- 
tion in the LTI phase has a sizable contribution from 
relativistic mechanisms 31, 32|. Furthermore, recent op- 



tical measurements revealed the presence of three elec- 
tromagnon peaks in the optical absorption spectrum of 
the non-collinear LTI state, the nature of which was not 
fully understood 



2 




h cj (meV) 

FIG. 1: Constant- Q scan on single crystals of YMn 2 Os at 
Q = (1,0,0) + qic at 4 K. Blue symbols are the neutron 
scattering data, while the red line is the result of spin-wave 
calculations. The inset is a contour map obtained from a 
powder sample of YM^Os at 1.5 K. 

In this Letter we report the results of inelastic neutron 
scattering (INS) for the LTI state of YMn 2 5 , provid- 
ing the first information about magnetic excitations in 
^LM^Os of multiferroic materials. We also present a 
theoretical analysis of the observed low-energy magnon 
bands, identifying them as electromagnon modes. Good 
agreement with results of optical experiments [l| clarifies 
the microscopic mechanism that couples spin waves to 
electric field in ^M^Os compounds. 

In order to study dynamic properties of YM^Os, a 
6 g powder sample and two single crystals with a total 
mass of about 1 g were used. The powder sample was sin- 
tered at 1180°C for 60 h with intermediate grindings. X- 
ray powder diffraction results at room temperature (not 
shown) indicated sin gle-p hase YMn 2 05 with a Pbam or- 
thorhombic unit cell [25j . Single crystals were grown us- 
ing E>203-PbO-PbF2 flux in a Pt crucible. The mixture 
was held at 1280°C for 15 hours then cooled to 950°C at 
1°C per hour. Both crystals used were characterized by 
neutron diffraction as described in Ref. |3lj |. INS experi- 
ments on the powder sample were carried out in a liquid 
4 He cryostat using cold neutrons with the incident energy 
of 25.25 meV on the Disk-Chopper-Spectrometer (DCS) 
at the NIST Center for Neutron Research (NCNR). INS 
experiments on the single crystals were performed on 
the cold neutron triple-axis spectrometer, SPINS, and 
on the thermal neutron triple-axis spectrometer, BT7, at 
NCNR. The single crystals were co-aligned in the (h, 0, 1) 
scattering plane within 1° and cooled in a 4 K closed- 
cycle refrigerator. For the INS measurements at SPINS, 
the energy of the scattered neutrons was set to Ef=5 
meV, and a liquid nitrogen cooled Be filter was placed be- 
fore the analyzer to eliminate higher order neutron con- 
tamination, while for the measurement at BT7, Ef=14.7 
meV and a PG filter were used. 

The inset of FIG. [T] shows a contour map of the INS in- 
tensity as a function of Q and Tiuo, obtained from the pow- 
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FIG. 2: Contour maps of the INS intensity, obtained at 
SPINS, along the (a) H- and (b) L-directions at 4 K. The in- 
tensities in (a) and (b) contours were multiplied by Vhui to en- 
hance the upper band. Corresponding calculated dispersions 
and intensities for parameters shown in Table [J convoluted 
with the instrument resolution ((c) and (d)). Experimental 
data taken at BT7 (open circles) and calculated intensities 
(contour map) for higher energy excitations from 10 to 20 
meV along the (e) H- and (f) L-directions. The radiuses of 
the circles are proportional to the integrated intensities of the 
observed peaks. 



der sample at 1.5 K. In the LTI phase, there exist three 
broad bands of excitations: one is centered at ~ 4 meV 
with a bandwidth of ~ 5 meV, indicating that it con- 
tains more than one excitation mode, and the other two 
modes, which are weakly dispersive, are centered at ~ 11 
and 16 meV. FIG. [U shows typical energy (ftoj) scans at 
Q = (1, 0, 0) + qic in the LTI phase. Up to ?iuj = 7 meV, 
there are at least three spin wave modes, resulting from 
a fit to centered at 0.98, 2.71, and 4.75 meV. We have 
performed similar Tiuo scans at numerous wave vectors to 
map out the dispersions of the spin waves along two high 
symmetry directions, (#,0,0.288) and (1.48,0,L). The 
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-0.2 
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0.15 


3.75 


2.75 


0.72 


0.13 


0.11 



TABLE I: The parameters of the spin Hamiltonian Eq.(pQ). 



contour maps in FIG. [2] (a) and (b) summarize our re- 
sults up to huj = 7 meV. Further constant- Q scans were 
performed at high energies around Tiuj = 15 meV. The 
data were fit by Gaussians and the peak positions are 
plotted as open circles in FIG.[2](e) and (f). As expected 
from the powder data shown in the inset of FIG. [TJ for 
Tiuj < 8 meV there are a few highly dispersive excitation 
modes, while for Tiuj > 8 meV the excitations are weakly 
dispersive. In addition, there is a gap of ~ 1 meV at the 
magnetic Bragg reflection due to magnetic anisotropy. 
The magnetic excitations are more dispersive along the 
i^-direction than along the L-direction, which indicates 
that the magnetic interactions are stronger within the 
afr-plane than along the c-axis. 

Since the spin structure of the LTI phase is rather com- 
plex, we approximated the magnetic structure to a flat 
spiral in the ab plane. We have numerically found the 
minimal-energy spin configuration in the magnetic unit 
cell corresponding to 39 x 1 x 4 crystallographic unit cells 
(i.e., replacing q IC = (0.48,0,0.288) by q' IC = (||,0, \)) 
using the Heisenberg spin Hamiltonian for a spiral state 
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where the sum is over pairs of Mn spins and Jij are the 
exchange coupling constants. We consider five nearest- 
neighbor interactions, J\ to J5 (defined as in Ref. l25f ) 
as well as the next-nearest-neighbor coupling along the 
c-axis, Jq (see FIG. [3j). We also include the easy plane 
anisotropics, V { = Ai (A 2 ) on Mn 3+ (Mn 4+ ) sites. The 
optimal values of the coupling constants obtained by fit- 
ting the INS data are listed in Table HI The red lines in 
FIG. [U and color maps of FIG. [2] (c)-(f) show the cal- 
culated INS spectrum. Despite the approximation made 
for the magnetic structure, our calculations reproduce 
the basic features of the observed magnetic excitations. 
The agreement between the calculations and the data is 
better for the independence than for the Independence, 
probably due to the fact that q' IC and qic are closer to 
each other in H than L. 

The basic features of INS results with the wave vector 
Q produced by a spiral state can be understood by con- 
sidering a circular spiral in the co-rotating-spin frame 
which has one coordinate axis parallel to the average 
spin vector and another axis perpendicular to the spiral 
plane [33j . In this co-rotating spin- frame, the momentum 
transfer vector Q' is conserved. The magnetic excitations 
at Q' by the magnetic field normal to the spiral plane can 
be observed at Q = Q' in the lab frame, while the same 
excitations by an in-plane magnetic field can be observed 




FIG. 3: Crystal and magnetic structures in the LTI phase of 
YMn205 projected onto the (a) ab- and (b) frc-plane. Thin 
white lines indicate the antiferromagnetic chains along the a 
axis. Also shown are the exchange coupling constants Ji — Jq 
used to calculate the spin- wave spectrum. Black arrows indi- 
cate the spin directions for the numerically obtained ground 
state of the spin Hamiltonian Eq.(l). The light blue arrows 
in (a) represent the directions of the oscillating Mn spins in 
the optical phason mode. Green (red) ellipses show that the 
angle between neighboring spins is increased (decreased) by 
this mode. 



at Q = Q' ± qic • 

Using Eq. (pQ) and the co-rotating-spin frame, we found 
that there are in total six modes that give rise to low- 
energy magnon bands in the LTI phase of YM^Os. 
Three of them - a phason and two rotations of the spiral 
plane - are generic for spiral magnets. The phason or 
sliding mode is a Goldstone mode of any incommensu- 
rate spiral state, having zero energy irrespective of the 
magnitude of higher harmonics in the spiral. It involves 
rotations of spins in the spiral plane and is excited ei- 
ther by an out-of-plane magnetic field with Q = 0, or by 
an in-plane field with Q = ±qic- The other two low- 
energy spiral modes are rotations of the ab spiral plane 
around the a and b axes, and have been found in other 
spiral magnets [34[. They can be excited either by the 
corresponding in-plane magnetic fields, H a and #5, with 
Q = (and, hence, can be observed as AFM resonances 
in optical experiments) or by an out-of-plane magnetic 
field with Q = ±qic- In the absence of an in-plane mag- 
netic anisotropy, these two modes are degenerate. The 
in-plane magnetic anisotropy induces higher spiral har- 
monics, lifting the degeneracy. 

These three modes have 'optical' counterparts, which 
is a specific property of ^M^Os compounds due to the 
fact that spins in a crystallographic unit cell belong to 
two AFM chains along the a axis. In the optical pha- 
son mode, spins in neighboring AFM chains rotate in 
opposite directions (see FIG. [3] (a)). The peculiar lat- 
tice geometry of ^M^Os leads to an almost complete 
cancelation of exchange interactions between neighbor- 



ing AFM chains [26 



28| , resulting in near degeneracy of 
the 'optical' and 'acoustic' magnon branches. 

FIG. |U (a) shows the calculated INS energies and in- 
tensities of these six modes at Q = (^|, 0, \). At this 
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FIG. 4: Calculated INS intensities of the first three acoustic 
modes (blue bars): the phason (P), rotations around the a 
and b axes (R a and ife), and their optical counterparts (O, 
red bars) at (a) Q = (f ,0, \) and (b) Q = (1 + §,0, ±). 
Some modes which have almost zero intensities are not shown 
here. 



Q position, the INS intensities corresponding to 'acous- 
tic' modes (blue) are large, while the intensities corre- 
sponding to the optical modes are hardly visible. At 
Q = (1 + H, 0, |) the situation is the opposite (see FIG.|4] 
(b)), which is related to the ~ | shift between neigh- 
boring AFM chains. Comparing these results with the 
experimental data at Q = (1.48,0,0.288), we have iden- 
tified the hu ~ lmeV peak as the optical phason which 
corresponds to the peak at 7.2 cm -1 observed in optical 
absorption [l[. The two optical rotation modes account 
for the broad peak around 2.5 meV in our INS data. 
This interpretation is supported by the observation of 
two close peaks at 16 cm -1 and 21 cm -1 in optical ex- 
periments [l|. These three peaks observed in the optical 
experiment are electromagnons as they are excited by 
electric field of light, E\\b [lj. It was previously specu- 
lated that the Heisenberg exchange mechanism that in- 
duces the electric polarization P\\b in the CM pha se also 

The 



28, 2c 



couples the optical phason mode to E\\b 
corresponding invariant for the Pbam space group has 
the form 

HME = —gEb I(Sl,n - S2, n -a) " (Se,n + ^7,n) 

n 

+ (S3,n — ^4,n-b) ' (&5,n + Ss,n)} 5 (2) 

where the indices 1-8 label Mn ions in the unit cell (see 
FIG. [3]), while n labels unit cells. The electric field 
of light, E\\b, modulates interchain interactions along 
the 6-direction, e.g., Js^g —> Js^g + A J, while J4,8+6 — >■ 
^4,8+6 — A J, which induces the relative rotation of spins 
in neighboring AFM chains (see FIG. [3]). To prove this 
scenario we find the coupling constant g using the value of 
polarization in the CM phase, P& = 1000 fiC- m -2 , and 
calculate the spectral weight, 5, of the optical phason. 
This gives S = 194 cm -2 , which is close to the exper- 
imental value of 170 cm -2 [28[. For large deviations of 



the LTI spin structure from the flat spiral, Eq.(j2j) can 
also explain the excitation of the spiral rotation modes. 

In summary, we identified six low-energy magnon 
modes and determined the exchange interactions in 
YMn205. Three of the six magnon modes are common 
to all spiral magnets, while the remaining three result 
from magnetic frustration, which is the source of the 
complex magnetic behavior and ferroelectricity found in 
AM1CI2O5 compounds. The latter three modes were iden- 
tified with the electromagnons observed in optical con- 
ductivity measurements. 
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